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Thermal gravimetric analysis (TGA) was applied to investigate the co-drying kinetics of biomass (Bermuda 
grass and red pine) with lignite. The drying rate of the lignite was much higher than that of the pure 
Bermuda grass and red pine. The drying of samples presented two distinctive falling rate periods. The 
activation energy of the pure Bermuda grass and red pine were higher than those of the grass/lignite 
and pine/lignite blends, respectively. In the first falling rate period, the dominant mechanism functions 
for the lignite, red pine and their blends were 1 - (1 - a) 1/3 , and the functions were 1 -(1 - ct) 0 - 5 in the 
Bermuda grass and its blends, while in the second falling rate period, the functions of the blends and their 
parent samples were -ln(l - a). The drying kinetics of blends had no linear relation with blend ratios. 
The blending of the lignite with the grass or pine could improve the drying of the individuals. 

© 2013 Elsevier B.V. All rights reserved. 


1. Introduction 

Lignite has high moisture content of 25-60%, which consid¬ 
erably limits its applications. For example, when low rank coals 
are burnt in boilers, about 7-10% of energy is used to evaporate 
the moisture in fuel [1 ]. It was calculated that the optimization of 
the drying process in future lignite power plants may lead to an 
efficiency increase of 4-6% [2-4], The development of an efficient 
lignite drying process is a necessary step towards the implementa¬ 
tion of the existing lignite-fuelled power plants [5,6], Meanwhile, 
biomass is a renewable energy source that will provide environ¬ 
mental and economic benefits [7-10], Due to the high volatile 
matter and moisture content, burning biomass alone cannot give 
very high energy output compared with coal combustion [11-13], 
Co-firing biomass with coal in existing boilers to generate elec¬ 
tricity has been proposed as a near-term, low-cost way to use 
biomass for reducing C0 2 emissions, furthermore, the existing coal- 
fuelled power plants can also be used with very few modifications 
[13-16], However, biomasses with high inherent moisture con¬ 
tent are similar to the lignite, and they should be dried before 
co-burning [17,18], TG method has many advantages in drying, 
such as ease of operation, online recording of the drying process, 
and minimal material requirement. Chen et al. [19] obtained the 
drying characteristics of cotton stalk using the TG method by apply¬ 
ing isothermal procedures. Tahmasebi et al. [20] gained the drying 
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kinetics of an Indonesian low rank coal and found that the dry¬ 
ing kinetics was best represented by two stages, the constant rate 
stage followed by the rate decay stage, and that the water removal 
rate was dependent mainly on drying temperature and coal sam¬ 
ple size. Tahmasebi et al. [21] revealed the chemical structure 
changes of Chinese lignite during fluidized-bed drying in nitro¬ 
gen and air. Kim et al. [22] examined the drying behaviour of 
low-rank Loy Yang coal in a fluidized bed dryer, and found that 
the drying rate increased with temperature and fluidization veloc¬ 
ity increasing and relative humidity decreasing. Pakowski et al. 
[23] experimentally determined the effective diffusivity of water 
in brown coal in superheated steam at 200 °C and atmospheric 
pressure using slightly compressed pellets of cylindrical shape. 
Many studies mainly focused on drying of foods or grains below 
100°C in order to preserve and store [24-26], For high moisture 
content matrixes, liquid diffusion is supposed to be the main mech¬ 
anism of moisture transport in the first falling rate stage of drying 
while vapor diffusion is supposed to dominate in the second falling 
rate stage [27], There were many studies primarily focused on 
co-combustion or pyrolysis of coal and biomass [28-32], Many 
researchers have revealed some results on the drying kinetics of 
pure biomass and lignite, but little information is available to date 
on the potential impact of co-drying kinetics for biomass and lignite 
at medium temperatures, especially co-drying kinetics of lignite 
and Bermuda grass or red pine are not found in the literatures. 

In the present work, thermal gravimetric analysis (TGA) was 
applied to investigate the isothermal co-drying kinetics of biomass 
(Bermuda grass and red pine) and lignite. The effects of tempera¬ 
tures, sample kinds and blending ratios on the kinetics parameters 
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Nomenclature 

a weight loss rate 

A the pre-exponential or frequency factor (s -1 ) 

DTG drying rate (mg min -1 ) 

E a apparent activation energy (J mol -1 ) 

k the reaction rate coefficient 

M mass fraction of sample 

R the universal gas constant [J/(mol I<)] 

t time (s) 

T drying temperature (K) 

Subscripts and superscripts 
0 initial 

oo final 


and mechanism were also investigated. This research will provide 
fundamental data for understanding the co-drying kinetics of 
biomass with lignite, and scientific reference for the design and 
improvement of the boiler. 

2. Experimental facility and test samples 

A thermogravimetric analyzer (TGA/SDTA851, Mettler Toledo) 
with a precision of 0.001 mg was used to conduct the drying experi¬ 
ments. Pulverized lignite, Bermuda grass and red pine particle were 
selected, and the biomass samples were crushed into 0.2-0. 6 mm 
in a high speed grinder, then, an electronic scale (FA1004, Heng 
Ping) with a precision of 0.1 mg was used to measure the mass 
of each sample. The mixtures were obtained from a platinum cru¬ 
cible in which the lignite/Bermuda grass and lignite/red pine blends 
were homogenized with biomass weight percentages of 20,40,60 
and 80 wt.%, and the initial mass of each sample was maintained at 
10 mg ± 0.5 mg. A crucible with the test samples was placed inside 
the drying chamber, and the drying processes were carried out 
under nitrogen with 40 mL/min flow. The thermal analyzer heated 
the sample from room temperature (20 °C) rapidly to the set tem¬ 
perature (50, 70,90,110,130,150 and 170°C) and then stabilized. 
Proximate analyses of the parent samples were performed in the 
analyzer, detailed measuring method from the references [33,34], 
The proximate analysis was listed in Table 1 . The heating values in 
Table 1 were estimated by using the method in reference [35], The 
ultimate analysis results for the parent samples [36-38] were also 
listed in Table 1. 


3. Kinetics analytical method 

The dependence of drying rate on temperature can be deter¬ 
mined by the Arrhenius expression [39]: 

**=A-e-&f(a) (1) 

where E a LJ/mol] is the apparent activation energy required for the 
moisture migration; A [min -1 ] the frequency factor; t [min] the 
time; R [J/(mol K)] the universal gas constant; and T[K] the absolute 
temperature; a the weight loss rate of samples at any time and it 
can be calculated using the following equation 


a = 


m 0 -m 


(2) 


where m is the weight of samples at any time, mo the initial 
weight of samples, and the weight of samples at the end of the 


drying process, respectively. By the way, the mass fraction of sam¬ 
ples, longitudinal coordinate in TG curves can be expressed as 



The integration function of Eq. (1) is given by: 

g{cl)= fm (4) 

Substituting of Eq. (1) into (4), and based on the isothermal 
kinetics, the integration function is described as follows: 

g{a) = kt+b (5) 


where k is the reaction rate coefficient, g(a) a function depending on 
the mechanism assumed for the reaction. Substituting each com¬ 
mon mechanism function shown in Table 2 into Eq. (5 ), and then the 
function is obtained by using the double logarithm method based 
on Avrami-Erofeev equation [40-42]: 

1 - a = (6) 


Double logarithm form about Eq. (6) is written as: 


ln[-ln(l - a)] = min t-i- In k 


(7) 


By plotting ln[—In(l — a)] vs.ln t, the m value is obtained at 
certain drying temperature. The m value is compared with those for 
the different kinetic models in Table 2, then the only mechanism 
function is selected. Furthermore, k value is got by plotting g(a) vs. 


The logarithm form of Arrhenius principle is set down [43]: 

lnk = -|| + lnA (8) 

By plotting In (k) vs. 1/T, E a and A can be subsequently related 
to drying conditions. 


4. Results and discussion 

4.1. Thermogravimetric and differential thermogravimetric 
analysis 

The TG and DTG curves of samples in different temperatures 
were obtained, but only partial results were presented due to the 
limited space. Figs. 1 and 2 illustrated the TG and DTG curves, 
respectively. The drying rate was equal to DTG value which was the 
variation of sample weight per unit time during the drying process. 

From the figures, the drying rate of the lignite was much higher 
than that of the pure Bermuda grass and red pine at the same 
drying temperature, and those of the Bermuda grass/lignite and 
pine/lignite blends were between the lignite and Bermuda grass 
or red pine, respectively. The drying of the Bermuda grass and red 
pine presented thermal lag characters compared with that of the 
lignite. This could be due to much bound water (bound to the inner 
pores of raw materials) in biomass and much free water (bound to 
the particle surface and existing in pores of raw materials) in the 
lignite with strong water-absorbing. The drying rates of the pure 
Bermuda grass and grass/lignite blends were higher than those 
of the pure red pine and pine/lignite blends pure due to that the 
red pine was belong to woody biomass with dense structure, and 
low moisture contents as shown in Table 1. As illustrated in Fig. 2, 
drying of samples presented three distinct periods, initial period, 
first falling rate period and second falling rate period. The initial 
period was the beginning of drying, which drying rate depended 
on the external process parameters. The slope of the DTG curve 
increased in a short time and the drying rate reached its maximum 
value rapidly. In this period, evaporation occurred only from the 
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Table 1 

Proximate and ultimate analyses of samples. 


Samples Proximate analyses (%) 


Ultimate analyses (%) 


HHV (kj/kg) 


M ad a V ad FC ad Art C„ b H„ JV ar O ar S ar 


Lignite 9.73 30.96 47.39 11.92 56.72 3.61 0.71 15.98 9.73 21.67 

Bermuda grass 7.26 67.82 11.52 13.4 46.80 6.12 0.91 41.66 0.25 17.91 

Red pine 5.35 64.66 21.65 8.14 49.30 6.30 0.01 44.10 <0.01 17.67 


a ad refers to moisture content, volatile substance content, fixed carbon content and ash content in sample corresponding to air dried basis. 


Table 2 

The common mechanism functions for solid reactions. 


Functions g(a) m Functions 


Di(a) <x 2 = kt 

D 2 (a) (1 - a)ln(l - a) + a = kt 

D 3 (a) p-(i-«) ,/s f = kt 

FiW— ln(l —a) = kt 


RiM 

RsCa) 


M«) 

Met) 


g(ct) 


1(1 -ct) V2 = kt 
1-(1 -a)' /3 = kt 
[—ln(l -a)] 1 ^ 2 = kt 
[— ln(l -a)] 1 ^ 3 = kt 


surface of the sample and since heat transfer promoted the evap¬ 
oration, the radius of the liquid/vapor interface over the exterior 
surface decreased continuously. Meanwhile, the rapid heating may 
have caused considerable expansion (puffing) of the small sam¬ 
ples, resulting in rapid evaporation and vapor diffusion. This period 
existed for a very short period of time, approximately 0.5-1.5 min 
depending on the drying temperature and blending ratio used.The 
maximum drying rate point was taken as a turning point (the first 
critical point) of the initial period and the falling rate period [ 17,43 ]. 
The first falling rate period, just beyond the first critical point, the 



(a) Bermuda grass / lignite at 130"C drying temperature 



(c) 80G20L a 


liquid/vapor interface had entered the pores and the liquid receded 
into the pores. During this period, the rate of evaporation of sur¬ 
face of outside decreased rapidly due to the increasing resistance 
of interface moving towards internal and the decreasing moisture 
content, the drying rate clearly decreased. During the transition 
of the first falling rate period to second falling rate period, there 
was possibly a change from porous tissue structure to an increas¬ 
ingly non-permeable structure. The shrinkage might have caused a 
denser and tougher structure which might increase the resistance 
to internal moisture movement, therefore the drying rates in the 



(b) Red pine/lignite at 130'C drying temperature 



(d) 80P20L b 


Fig. 1. TG i 


; of biomass/lignite. 
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(a) Bermuda grass /lignite at 130"C drying temperature (b) Red pine / lignite at 130'C drying temperature 






(c) 80G20L ( d ) 80P20L 


Fig. 2. DTG curves of biomass/lignite. 


second falling rate period decreased slowly than the first falling rate 
period. Casehardening was another possible reason for low drying 
rate in the 2nd falling rate period. 

The drying rate of the blends increased with increasing the dry¬ 
ing temperature because of the reinforce of moisture migration, 
which resulted from the enhancement of driving force on heat 
transfer in the samples. For example, the maximums of the drying 
rates at 170°C were about 16 times and 15 times higher than those 
at 50 °Cfor the 80G/20L blend and 80P/20L blend, respectively. 


4.2. Kinetic analysis 

In order to accurately reveal the co-drying behaviour of blends 
in whole falling rate stage, the drying kinetics of the first falling 
rate period and the second falling rate period were distinguished. 
Fig. 3 presented the linear fitting by plotting In [- In (1 - a)]~ In t 
for the 80G20L blend in the first falling rate period at 130 °C. The 
m value was 1.091, which was close to 1.11, so f?2(«) was taken 
as a suitable mechanism function. The k value was 0.394 min -1 , 
which was obtained by plotting the g(a) vs. t, as shown in Fig. 4. 
The method might also be applied to other drying temperatures, 
and other period, and other samples. 

The activation energy and the frequency factor of blend were 
obtained by plotting In (k) vs. 1 /T in the first falling period as shown 
in Fig. 5, from which E a and A values were predicted, 15.08 kj/mol 
and 10 1 - 53 min -1 , respectively. The method might be also applied 


to other period and other samples. The results were presented in 
Table 3. 

The activation energy of the lignite was higher than those of the 
pure Bermuda grass, red pine and blends. The activation energy 
of the pure Bermuda grass and red pine were higher than those 
of the grass/lignite and pine/lignite blends, respectively, indicat¬ 
ing that the blending of the lignite with the Bermuda grass or red 



lnt 

Fig. 3. Linear fitting of drying of the 80G20L blend at 130 °C by lnln method. 
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Time (min) Ratio (%) 


Fig. 4. Linear fitting between G(a) vs. t for the 80G20L blend at 130 °C. 


(a) Trailing rate period 



l/T 

Fig. 5. Arrhenius-type relationship between In (k) vs. 1 /T for the 80G20L blend in 
the first falling period. 



Ratio (%) 

(b) 2 nd falling rate period 

Fig. 6. Comparison between calculated activation energy and real values. 


pine could improve the drying of the lignite and Bermuda grass 
or red pine. The grass/lignite and pine/lignite blends can be dried 
under lower temperature than their parent samples. The activation 
energy of samples in the first falling rate period was much higher 
than that in the second falling rate period, which was attributed 
to the reducing of bound water in the sample. In order to find 
out whether interactions between the components of the blends 
occurred, the activation energy of the blends calculated by using 
the ratio of samples in the first falling rate period and the second 


falling rate period as shown in Fig. 6. The drying kinetics character¬ 
istics of the blends had no linear relation with blend ratios due to 
interaction between the lignite and Bermuda grass or red pine, and 
the real values were much lower than the calculated values, indi¬ 
cating that there were synergistic interactions between the parent 
samples during the co-drying of the lignite and Bermuda grass or 
red pine. 


Table 3 

Kinetic parameters and mechanism functions of the drying of samples. 


Samples 


1st falling rate 


100L 

20G80L 

40G60L 

60G40L 

80G20L 

100G 

20P80L 

40P60L 

60P40L 

80P20L 

100P 


E a (kj/mol) A (mim 1 ) g(a) R SD E a (kj/mol) Afmin- 1 ) 


18.53 
15.85 

14.53 
13.87 
15.08 

16.54 
16.51 
14.79 
14.35 
15.45 
17.22 


10 1 - 82 R3 -0.999 

10 1 - 61 R2 -0.994 

10 1 - 52 R2 -0.995 

10 1 - 46 R2 -0.997 

10 1 - 53 R2 -0.999 

10 1 - 74 R2 -0.999 

10 1 - 61 R3 -0.997 

10 1 - 36 R3 -0.995 

10 1 - 30 R3 -0.999 

10 1 - 53 R3 -0.996 

10 1 - 79 R3 -0.999 


0.006 15.26 10 1 - 74 

0.011 12.96 10 1 - 49 

0.008 11.25 10 1 - 28 

0.004 10.67 10 1 - 21 

0.008 12.27 10 1 - 33 

0.003 14.1 10 1 - 67 

0.005 13.31 10 1 - 53 

0.004 12.19 10 1 - 36 

0.006 11.96 10 1 - 28 

0.008 12.7 10 1 - 45 

0.006 14.93 10 1 - 68 


Sfce) 


-0.998 

-0.998 

-0.997 

-0.997 

-0.998 

-0.997 

-0.996 

-0.997 

-0.996 

-0.999 

-0.997 


SD 

0.006 

0.003 

0.005 

0.004 

0.004 

0.006 

0.008 

0.006 

0.007 

0.008 

0.005 
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Moreover, in the first falling rate period, the mechanism func¬ 
tions of the drying for the lignite, pure red pine and their blends 
were 1 - (1 - cr) 1 ^ 3 , which was controlled by the phase boundary 
reaction and spherical symmetry. The rate of water loss is domi¬ 
nated by the rate of an interface advance, which generated between 
reactant and product phase. The formula is also known as the con¬ 
tracting volume (cube or sphere) equation, which can describe rate 
variation of interface advance in different directions. The mecha¬ 
nism functions of the Bermuda grass and its blends with the lignite 
were 1 —(1 — a) 0 - 5 , which was controlled by the phase boundary 
reaction and cylindrical symmetry. This rate equation derived from 
a reaction readily initiated at all surfaces of a cylinder or column¬ 
shaped prism of reactant, or the interface moving inwards from 
the edges of a disc or plate-like particle. In the second falling rate 
period, the functions of the grass/lignite blends, pine/lignite blends 
and their parent samples were -ln(l — a), which was controlled 
by the random nuclear producing and growing process. Nuclea- 
tion is often considered as a series of chemical changes across 
reactant surfaces, which initiated by an advancing energetic wave 
occurring randomly at the reactant/product interface. Growth is the 
maintained reaction within the interface between the reactant and 
product phases. 

5. Conclusions 

The drying rate of the lignite was much higher than those of 
the pure Bermuda grass and red pine, and those of the grass/lignite 
and pine/lignite blends were in between those of the lignite and 
Bermuda grass or red pine. The drying of the Bermuda grass and 
red pine presented thermal lag characters compared with that of 
the lignite. The drying rate of blend increased with increasing the 
drying temperature. The drying of samples presented two obvious 
falling rate periods. In the first falling rate period, the mechanism 
function of drying for the lignite, pure red pine and their blends 
were fitted to 1 - (1 - a) 1 ^ 3 , and the functions ofthe Bermuda grass 
and its blends with the lignite were fitted to 1 -(1 - a) 0 - 5 , while 
in the second falling rate period, the functions of the grass/lignite 
blends, pine/lignite blends and their parent samples were fitted to 
-ln(l -a ).There were synergistic interactions between the parent 
samples during the co-drying of the lignite and Bermuda grass or 
red pine. 
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